Treatment of Correlation Effects in Electron Momentum Density: 
Density Fuctional Theory and Beyond 



B. Barbiellini and A. Bansil 

Department of Physics, Northeastern University, Boston, MA 02115 USA 

Abstract 

Recent high resolution Compton scattering experiments clearly reveal that 
there are fundamental limitations to the conventional local density approx- 
imation (LDA) based description of the ground state electron momentum 
density (EMD) in solids. In order to go beyond the framework of the density 
functional theory (DFT), we consider for the correlated system a BCS-like 
approach in which we start with a singlet pair wavefunction or a 'geminal' 
from which the many body wavefunction is then constructed by taking an 
antisymmetrized geminal product (AGP). A relatively simple practical im- 
plementation of the AGP method is developed where the one-particle orbitals 
are approximated by the Kohn-Sham solutions used in standard band compu- 
tations, and the orbital-dependent BCS energy scale Aj is determined through 
a readily computed exchange-type integral. The methodology is illustrated by 
considering EMD and Compton profiles in Li, Be and Al. It is found that in 
Li the present scheme predicts a substantial renormalization of the LDA re- 
sult for the EMD; in Be, the computed correlation effect is anisotropic, while 
in Al, the deviations from the LDA are relatively small. These theoretical 
predictions are in qualitative accord with the corresponding experimental ob- 
servations on Li, Be and Al, and indicate the potential of the AGP method 
for describing correlation effects on the EMD in wide classes of materials. 
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I. INTRODUCTION 



Recent high resolution Compton scattering experiments on a number of materials clearly 
show that the EMD of the ground state is not described satisfactorily by the conventional 
LDA-based framework In Li, the size of the discontinuity in the EMD at the Fermi 



momentum -pp appears to be anomalously small Q -nearly zero 0-^, and differs sharply 
from the values of 0.55 — 0.85 deduced via a variety of studies of the correlated homoge- 
neous electron gas (HEG) stretching over last several decades [|T^T7|. In Be, comparisons 



between highly accurate computed and measured Compton spectra |0] indicate systematic 
discrepancies which are anisotropic and would be difficult to explain within the standard 
Lam-Platzman (LP) type correlation correction to the EMD |1T8|-pO|l which is isotropic by 
construction. The EMD in Al, on the other hand, is adduced to be reasonably close to the 



LDA predictions [p|-[TO||. Several attempts have been made to gain a theoretical handle on 



this problem pT|-p^ . The early optimism of a GW computation in Li did not hold up to 



a later investigation p2| , p3| . Surprisingly, Quantum Monte-Carlo (QMC) studies of Si and 
Li p^ , |25[| also do not reveal any substantial differences in the EMD with respect to the LDA 
predictions. It should be emphasized that a successful theory must not only explain what 
may be thought of as the "low" density case of monovalent Li, but also it must approach 
the LDA results correctly in the "high" density limit of trivalent Al in addition to explain 
the behavior of divalent Be. 

In order to make progress, we observe first that the current GW and QMC computations 
in connection with the EMD assume that the fermionic correlation is not modified drastically 
by interactions and that an adiabtic path exists between the free and the interacting electron 
gas. The many body wavefunction underlying the GW and QMC work is built from Slater 
determinants of single-particle orbitals, with an implicit nodal structure and plasmon-type 



^We will not concern ourselves here with possible complications in the interpretation of Compton 
data with regard to questions of background subtraction, failure of impulse approximation, etc. 
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physics of correlations p6[ which is more or less similar to that of the LDA [^. Bearing 
these considerations in mind, we have been motivated to examine other strategies, and here 



we propose a BCS-like approach for going beyond the framework of the DFT [p8|-p3|], wherein 



we start with singlet electron pairs- referred to as 'geminals' |^^, and then combine these 
objects into an antisymmetrized product (AGP) to obtain the wavefunction of the many 
body system. Such an idea has been invoked previously, for example, in connection with 
liquid He^ and in quantum chemistry calculations of molecules [^. An application to 
the problem of the EMD in extended systems which is presented here has, to our knowledge, 
not been reported in the literature before. Notably, the AGP wavefunction in molecules 



typically yields only 40 — 50 % of the available correlation energy |^0|]. However, our main 
goal here is to develop a working scheme which can capture some of the essential physics 
of the EMD of the correlated electron gas, and not so much to obtain the total energy 
accurately. In this spirit, it is hoped that the AGP method can provide a useful tool for 
understanding the recent Compton spectra and the nature of the EMD over a wide range 
of electron densities. 

An outline of this article is as follows. We start in Section II with a rigorous expression for 
the EMD, p(p), in the correlated electron gas in terms of the eigenvalues and eigenfunctions 
of the one-particle density matrix p(r, r') |]35| , |36| . Section III discusses the conventional 
Lam-Platzman (LP) treatment [[l^ -[20[| of correlations within the LDA which is based on the 



momentum density Ph{p) of the interacting HEG p^-[T7|, including a useful parametrization 
of the more recent QMC data for Ph{p) fl^l- Section IV takes up a description of the AGP 



wavefunction and the associated total energy functional A relatively simple practical 
implementation of the AGP scheme is discussed in Section V. The theory is illustrated by 
considering aspects of the EMD and Compton profiles in Li, Be and Al in Section VI. We 
emphasize that although these results indicate that the present method is promising, further 
work is necessary for confronting the AGP model quantitatively with experiments. Finally, 
Section VII summarizes the main results of this article. 
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II. MOMENTUM DENSITY AND NATURAL ORBITALS 



We start with the one-particle density matrix, /3(r,r'), which is defined in terms of the 
normahzed N-particle wavefunction, as 

p(r,r') = iv/cie^*(r,0^(r',0, (1) 

where the integral extends over the coordinates of all other particles. The eigenf unctions 
and eigenvalues of p define the natural orbitals, ipi, and the associated occupation numbers, 
rii [0. The density operator p can then be expanded into projectors \ ipi >< ipi \ via the 
spectral theorem: 

oo 

p = ^ni\'4)i>< ipi\ . (2) 

i=l 

If one requires the natural orbitals to possess the symmetry of the Hamiltonian, then they 
constitute a unique decomposition in an ortho normal basis set. Other basis sets for describ- 



ing the many body wavefunction are often used, e.g. the generalized overlap amplitudes p6 
which are linearly dependent. Any basis set can of course be related to the natural orbitals 
through an appropriate canonical transformation. 

The EMD is defined as the momentum transform of the density matrix. 



P(P) = 7 J drd r p(r, r) exp(-zp ■ (r - r')) . (3) 

Note that, in general, the EMD involves off-diagonal elements of the real space density 
matrix. In terms of the natural orbitals ipi, however, the EMD can be cast in the simple 
form 

P(P) = E'^i I < P I > (4) 

i 

where < p \ ipi > is the momentum transform of ipi. 

If the many body wavefunction is represented by a single determinant, which is true in 
the case of Hartree Fock or the DFT then the density matrix is idempotent (p = p^) 



and reduces to a summation over the occupied spin-dependent orbitals ipi, i.e. 



N 

p(r,r')=E^«(rM(r')- (5) 

i=l 

The EMD in this independent particle model (IPM) is similar to the more general Eq. (2), 
except that the occupation numbers are now strictly or 1 depending upon whether the 
state in question is empty or filled. Electron correlations, in effect, allow these occupied 
and empty IPM states to mix so that the occupation of states below the Fermi energy (Ep) 
becomes less than 1, while that of states above Ep becomes non-zero. 

One further point is noteworthy. In a periodic system, the electronic states are the Bloch 
states, ipkfj.lr) , where /i is the band index and k is the crystal momentum which is restricted 
to the first Brillouin zone (BZ). The associated occupation numbers, = N^{]s.), possess 
translational symmetry, i.e. 

iV^(k) = iV^(k+G„), (6) 
with respect to the set of reciprocal lattice vectors G„. 



III. CORRELATION CORRECTION IN THE DENSITY FUNCTIONAL THEORY 



As a consequence of the Hohenberg-Kohn theorem ||2^, the ground-state expectation 



value of any operator O is a functional 0[n(r)] of the electron density n{r). It can be shown 
that [g 

0[n{r)] = Oo[n{r)] + ^E^Mr)m , (7) 

where Oo[^(r)] is the expectation value for a system of noninteracting fermions moving in the 
effective field involved in the Kohn-Sham formalism, E^^ is the exchange-correlation energy 
functional and A is a scalar coupling parameter for the operator O. The value obtained 
for the non-interacting case- given by the first term on the right side of Eq. (^), must be 



corrected by the second term. This so-called LP correction [18 1 provides a formal scheme 



for treating correlation effects in the interacting system within the DPT. 
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Insofar as the momentum density in extended systems is concerned, LDA is perhaps the 
most common implementation of the DFT. The relevant expression for the LP correction, 



Ap[n(r)](p), is based on the treatment of E^c in the HEG [|T9 



Ap[n{r)]{p) = [ph{rs{r),p) - po(p)] ^^(r) d\ . (8) 

Here, the integral extends over the Wigner-Seitz (WS) unit cell. The square brackets in the 
integrand give the difference between the momentum densities ph and po of the interacting 
and non-interacting HEG, respectively, evaluated at the local density n(r) of the physical 
system given by the electron density parameter Vgij). 

Despite the limitations of form (H) -a point to which we return below, this equation has 
been used extensively in the literature for evaluating the LP correction by using a variety of 
different results for the HEG starting with the work of Daniel and Vosko [p!2| -[T7[|. The most 
recent QMC data for ph[j>) however would be the most satisfactory to employ [jl7 . 



With this motivation, we provide a convenient formula for pySjs^V) as a function of r^, 
which we have obtained by parameterization of the QMC data []T^. For HEG of density n 

We write QMC results for ph{v) as 

r a(l) (1 -a(2) a;2), if a; < 1; 

Ph{r,,v) = \ (10) 
[ a(3) exp(— a(4) [x — 1)) + T/x* otherwise, 

where x = p/pr with pp = (9/47r)^/^/rs, and a(l) through a(4) are r^-dependent fitting 
parameters. Our form of ph{rs,p) differs somewhat from that given by Farid et al. [jl6 



m 



that the tail in Eq. (|ToD is exponential while Ref. uses a Gaussian tail. The coefficient 
T is given by Yasuhara and Kawazoe 



as 



T = U^?9iO), (11) 



where a = (4/97r)^/^ and 
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64 

" (8 + 5 r, + 383/1200 rlf ' ^^^^ 
is the pair correlation function at r = 0, which has been obtained analytically by Overhauser 
PP| . Finally, we have fitted the coefficients a(l) through a(4) to the QMC data and obtained 
the following values. 

a(l) = 1 - .010 r„ 
a(2) = 0.025 r„ 
a(3) = 32/13 (5, 

a(4) = 4, (13) 

where 5 

5 = l/3-a(l) (l/3-l/5a(2))-^ (14) 



has been chosen to satisfy the electron number sum rule. Fig. |1| shows that the QMC results 
parametrized here over the typical metallic range of densities (r^ varying between 2 and 5) 
differ significantly from those given by the early work of Lundqvist ||T3[ parametrized by 
Cardwell and Cooper |2^; the original QMC data points lie indistinguishably close to the 



fitted solid curve of Fig. 1 and are not shown for simplicity. 



IV. TREATMENT OF CORRELATION EFFECTS BEYOND THE DFT 

The LP expression of Eq. (|^) suffers from a number of limitations. Most importantly, 
this form only describes an isotropic redistribution of the EMD, even though the effect of 
correlations will in general be anisotropic. This also implies that the LP corrected EMD will 
possess a spurious feature at the free electron Fermi radius in all directions in the momentum 
space. Finally, the movement of spectral weight from below to above Ep in Eq. (^ is not 
very transparent in the sense that it is not associated directly with changes in the occupation 
numbers of various states, making it very difficult to impose reasonable physical conditions 
on the occupation numbers such as those of translational symmetry of Eq. (|^). 
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Bearing these considerations in mind, and the fact that the QMC and GW results based 
on the LDA in Li p3| , |25| are not all that different from the LP correction and, therefore, 
would not explain the experimental EMD in Li, it is clear that there is need to go beyond 
the basic form of the Slater- Jastrow many-body wavefunction implicit in much of the QMC 
work given by 

^ = FD^D^, (15) 
where the Jastrow factor F accounts for the plasmon zero modes p6| and and are 



Slater determinants for up and down spin electrons. To make a headway, we propose using 
a BCS-like many body wavefunction in which individual elements involve singlet electron 
pairs, rather than one-particle orbitals. The underlying nodal structure of such a many body 
wavefunction tends to minimize the effects of the exclusion principle ||28|. As already noted. 



such an approach has been invoked in a variety of problems in the literature over the years 
P(]|^3[, and here we are suggesting its use in connection with the treatment of momentum 
density in extended systems. 

The starting point is a singlet pair wavefunction or the associated generating germinal 
0(ri t, Y2 i) |32,34|, which can be expressed in a diagonal expansion of natural spin orbitals 



with undetermined coefficients gi?,: 

^(ri T,r2 i)=Y. 9r ^*(ri)^.(r2) [| T> I i> "I i> I !>]• (16) 

i 

The many body wavefunction then is the AGP 

* = i{nSV(r.T,r,i)} , (17) 

where A is the antisymmetrization operator. The AGP belongs to the category of wavefunc- 
tions of extreme type Eq. (|1^) can be cast into an N/2 x N/2 determinant 



^ Incidentally, if N is odd, the extreme type wavefunction is proportional to the antisymmetrized 
product of a one body orbital / and {N — l)/2 factors of 4>] f must be strongly orthogonal to (j), 
that is, J f{ri)(j){ri,r2)d^ri = for all [||. 
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^ = Det|0(r, t,r, (18) 
The one-body wavefunctions i/ji are the natural orbitals for not only the generating geminal 



(f), but also for the AGP wavefunction ^ For a two electron system, AGP is exact and 
equivalent to a configuration interaction calculation. 

The ipi^s and Qi^s may be determined by minimizing the total energy functional. In 
extended systems the coefficients Qi are not particularly convenient, and it is more useful to 
introduce a new set of coefficients hi to define a " Cooper pair" function given by: 

C(ri T,r2 i) = E ^r(ri)^.(r2) [| T> | i> -| i> | T>] (19) 



The non-Hartree-Fock-like term in the total energy can then be written as ||3^ 



E[h,, V^,] = Ehf[p] + Ebcs[K, + 0{l/N) , (20) 



where Ehf is the Hartree Fock functional and Ebcs is a BCS-type functional [p9ip2| given 

by 

Ebcs = \<C\Vi2\C> . (21) 

The normalization of the AGP wavefunction imposes a relationship between the coefficients 
hi and the occupation numbers rii via the condition p2[ 



hi = ±^ni{l-ni). (22) 

Note that for Coulomb interaction the pair potential V12 is repulsive, so that energy can be 
gained only through the exchange part Ehf of the Hartree- Fock functional. It can be shown 
that the amplitudes hiS then must change sign at the "pseudo- Fermi surface" corresponding 
to the Hartree- Fock solution (i.e. the solution that neglects the term Ebcs) 113] ■ The energy 
can also be gained of course through the term Ebcs in Eq. ( pOf ) by the introduction of lattice 
dynamics as is the case at the superconducting transition pO[1 . 

Recently, Goedecker and Umrigar |^,^ have considered an approximation to the two- 
particle density matrix a which leads to a functional of natural orbitals similar to form 
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discussed here. However Csanyi and Arias have shown that the electronic states 



predicted by Ref. are overcorrelated; the reason is that a is varied over too large a class 
of functions without the restriction of A^-representability. This problem is circumvented in 
our case since the the AGP functional is iV-representable by construction. 

V. A SIMPLE IMPLEMENTATION OF THE AGP MODEL 

Our main goal is to gain a handle on the nature of occupation numbers in the corre- 
lated electron gas. In this spirit, we start by approximating the natural orbitals ipi by the 



Kohn-Sham orbitals for the one-particle crystal potential; the associated eigenvalues are 



denoted by Ei with Si = Q defining the Fermi level following the usual convention. 



The solution to the BCS problem involves an energy scale Aj which determines the 



mixing of states above and below the Fermi level. The occupation numbers are given by 

rii is obviously greater than 1/2 for Cj < 0, and less than 1/2 for Cj > 0; for e << Aj, nj — > 1 
and for e >> Aj, ni — > 0. At Si = 0, = 1/2. In order to make progress at this point we 
require a method for estimating the value of Aj. As already noted, there will be cancellation 
between the terms Ebcs and E^f in Eq- (pOD, and we would physically expect the total 
energy of the AGP wavefunction to be quite close to that of the Kohn-Sham Hamiltonian. 
With this motivation, we make the assumption that the energy cost of electron pairing will, 
in fact, he roughly compensated by the exchange energy. If so, it readily follows that 

|A,|~/, , (24) 



where 



J on ri on ro 



(25) 



' 6n{Yi)6n{r2) 

is an exchange type integral . We emphasize that the validity of assumptions leading to 
Eq. ( P^ ) is unclear and that this point deserves further investigation. Within the LDA, an 
explicit formula for Jj can be obtained straightforwardly: 
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where ^^^(r) = 2/7r(37r^r?,(r)^/^ is the Kohn-Sham exchange potential As expected, J, is 



large when tpi{r) is confined to a small spatial region. Finally, the n^'s are renormalized in 
order to satisfy the correct sum rule, i.e. 

iV = E^- (27) 

i 

The evaluation of the EMD follows along the lines of the standard IPM computations, except 



that the occupation numbers for the correlated electron gas given by Eq. (|23D are used. For 
the specific EMD and Compton profile results on Li, Be and Al reported in Section VI below, 
we have used orbitals ipi calculated using linear muffin-tin orbital (LMTO) band 
structure method for convenience, including corrections to the EMD due to the overlapping 



sphere geometry 



The essential parameter which controls the correlation effect in our theory is, 

a = A/w ~ li/w, (28) 



where w denotes the relevant valence electron band-width. By keeping Eqs. (^3|) and 
in mind, it is evident that when Jj is comparable to w, a ~ 1, and states deep in the Fermi 
sea are renormalized inducing significant shift of spectral weight from below to above £j = 0. 
On the other hand, when li << w, a is small and only states near the Fermi momentum 
are redistributed. 

Eq. (26) shows that when ipi's are the Bloch wavefunctions ■i/'^k, then the Jj's and hence 
the Aj's will be energy and k-dependent, and the resulting correlation correction to the EMD 
will in general be anisotropic. Concerning the break at the Fermi momentum, solutions for 
some purely repulsive model interactions indicate that A^(k) can oscillate as a function of 
energy, and become antisymmetric with respect Ei, going to zero at the Fermi level. |^ If 



so, the present scheme admits the presence of discontinuities at the Fermi momentum in the 
EMD. Notably, the AGP satisfies the translational-symmetry requirement of Eq. (|^). 
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VI. ILLUSTRATIVE RESULTS ON Li, Be AND Al 



In this section, we briefly discuss some of tlie salient features of the EMD predicted by 
the relatively simple theoretical model of the preceding subsection. Examples of Li, Be and 
Al are considered as a means of getting a handle on correlation effects for a range of electron 
concentrations- varying from one electron per atom in Li to two in Be and three in Al. 
As emphasized already in the introduction, our main purpose here is to establish that the 
present scheme is capable of explaining, at least in principle, some of the key experimental 
observations in Li, Be and Al P|-pD|]. 

In connection with the Compton spectra, it should be noted that the Compton profile 
(CP), J{pz), represents the double integral of the ground-state EMD p(p): 



where Pz lies along the scattering vector of the x-rays. 

Li is considered first. Using Eq. (p6D, we obtain ^ 0.1 Ry, which when compared with 
the valence electron band width, w ~ 0.25 Ry, yields a ^ li/w = 0.4. This implies that 
the IPM values of occupation numbers will be modified substantially by correlations. Fig. 
2 compares a typical CP (results along [100] are given) in Li for LDA with and without the 
LP correction, and the AGP model computation. The LDA (solid) displays a cusp around 
0.6 a.u. which reflects the break in the EMD at pp- The addition of the LP correction 
(light dashed) shifts some spectral weight from low to high momenta, but the cusp at pp 
is essentially unchanged. In sharp contrast, in the AGP curve (thick dashed), the feature 
at piT' is completely smoothed out and the EMD is redistributed throughout the momentum 
space. The AGP results are at least qualitatively similar to the behavior of the experimental 
CP's in Li which show that the size of the discontinuity Zk at pp is very small and that 
the measured CP is substantially lower than LDA around p = We should keep in 

mind nevertheless the fact that the AGP results in Fig. 2 invoke a series of approximations 
which might oversmooth the discontinuity at pp. Notably, spin-polarized electronic structure 




(29) 
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computations indicate the importance of spin- wave fluctuations in Li suggesting that 



some renormahzation of the occupation numbers in Li could also be produced by such a 
mechanism. 

Interestingly, an extensive study of Liioo-xMgx disordered alloys over the composition 
range < x < 40 indicates that the measured CP's come into better accord with the 
LDA predictions with increasing Mg concentration. This observation will find a natural 
explanation in our theoretical framework since we expect the width of the valence band to 
increase with increasing Mg content with little change in the exchange integral /», so that 
the parameter a (see, Eq. ^8]) will presumably decrease as will the correlation effect. 

Turning to Be, the exchange integral is computed to be, Jj ~ 0.1 Ry, but w ~ 0.8 
Ry, giving a ~ 0.13. Therefore the correlation effect on the EMD in Be is smaller than Li, 
consistent with experimental observations Fig. 3 focuses on the difference, A J-^'^^°''^ = 

J AGP — JldAi between the AGP and the LDA (including the LP correction) predictions of 
CP's along two different crystal directions. The LP correction (dashed) is seen to mimick 
the overall shape of A J-^'**^"'"^. However, in contrast to the LP correction which is isotropic, 
the residual correlation correction given by AJ^heory amsotropic. Along [100], AJ^^^°''y 
differs relatively little from the LP curve, but along [001] the excursions are much larger. 
The results of Fig. 3 are in qualitative accord with some of the characteristics of the 
experimental CP's in Be. By comparing the residuals of our Fig. 3 with the experimental 
residuals (defined as, AJ^^^^ = JExpt — Jlda, which are given in Fig. 4 of Ref. [0), it will 
be seen that AJ-^^^*, like AJ™^"*"^, deviates little from the LP curve along [100]; in sharp 
contrast, along [001] A J^^^* contains a remarkable pattern of undulations rather similar to 
that seen in AJ'^^'^"^^ in Fig. 3. Further study is needed to pin down the origin of specific 
structures in AJ'^^'^"'^^ in terms, for example, of the Fermi surface features. 

Finally, in the case of Al, we find Jj ~ 0.05 Ry, which as expected is small, reflecting the 
decrease of Jj with increasing electron density. On the other hand, the band width w ^ 1 
Ry is quite large, so that a ~ 0.05 is rather small. Therefore, the EMD in Al will not suffer 
much redistribution beyond the LDA via electron correlations. This again is consistent with 
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experimental observations which show that the EMD in Al is described reasonably well by 
the conventional LDA picture p|-p^. We illustrate this point further by Fig. 3 where the 
LDA and AGP momentum densities p{p) are compared along a typical direction (note these 
are 3D momentum densities and not the CP's). The AGP and LDA are seen to be quite 
close overall, some rounding of the Fermi break in AGP notwithstanding. 

VII. SUMMARY AND CONCLUSIONS 

We discuss the theoretical treatment of the EMD in the correlated electron gas with 
the purpose of developing an understanding of recent high resolution Compton scattering 
experiments on a variety of materials. EMD can be expressed rigorously in terms of the 
eigenfunctions and eigenvalues of the one-particle density matrix, which play the roles of 
effective orbitals ipi^s - the 'natural orbitals', and the associated occupation numbers rij's, 
respectively. Much of the existing EMD work in crystals is based on the use of Bloch 
wavefunctions computed within the LDA and the related IPM values of occupation numbers 
(i.e., 1 for filled and for empty states). In the DFT, the expectation value of any operator 
(excepting the electron density n{r)) must be corrected with respect to its IPM value in order 
to include the effect of correlations in the interacting system. The standard calculation of 
such an LP correction to the IPM momentum density involves the momentum density Ph{p) 
in the interacting HEG. Although many attempts to evaluate ph in the HEG have been made, 
the recent QMC results are perhaps the most accurate. In this connection, we provide a 
convenient parametrization of the QMC data for ph as a function of the electron density 
parameter r^; this form of ph should be used for computing the LP correction. Note that 
the LDA-based LP correction possesses fundamental limitations, e.g., it can only describe 
an isotropic redistribution of the EMD. 

In order to go beyond the framework of the DFT, we propose using a BCS-like wave- 
function in which the starting point is an electron singlet pair or a 'geminal', and the many 
body wavefunction is constructed as an AGP. Although such an approach has been invoked 
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previously in a variety of problems, we are not aware of its application to treat the EMD 
in extended systems. The AGP wavefunction tends to minimize the effects of the exclusion 
principle and possesses a fundamentally different nodal structure compared to the Slater- 
Jastrow type many body wavefunctions based on single-particle orbitals which are implicit 
in the LDA and QMC work. We begin to explore the nature of the EMD in the AGP 
scheme in this article by approximating the natural orbitals by the Kohn-Sham orbitals. By 
making the assumption that the energy gained through the BCS-like term will be roughly 
compensated by the exchange-correlation part of the Hartree-Fock term in the total energy 
functional, we show that the BCS energy scale, Aj ~ Jj, where /j is a readily computed 
exchange-type integral. The key parameter which controls the redistribution of states in our 
theory then is, a ~ /j/w, with w being the valence electron band- width. 

Finally, we consider the application of the AGP method to discuss the EMD and Compton 
spectra of Li, Be and Al as illustrative examples. In Li, we find a ~ 0.4, implying a 
substantial renormalization of the LDA states, in essential accord with experimental results 
which indicate the break in the EMD to be renormalized to a nearly zero value. In Be, 
a ~ 0.13, and the theory predicts anisotropic correlation effects in the EMD which show a 
remarkable resemblance to the anisotropy of measured Compton profiles. In Al, we estimate, 
a ~ 0.05, so that the correlations yield relatively little modification of the EMD, and 
here again, this is consistent with the experimental observations which indicate that Al is 
described reasonably within the conventional LDA picture. Taken together, these results- 
involving a range of electron concentrations, show the potential of the AGP method in 
providing a theoretical method for describing correlation effects on the EMD in wide classes 
of materials, although further work is necessary in this connection. Other open questions 
concern the excitation properties of the AGP ground state. 
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Figure Captions 



FIG. 1. Momentum density p(p) in the homogeneous electron gas. The momentum p is 
given in scaled units of the free electron Fermi momentum pp. Parameterized QMC (solid) 
and LDA (dashed) results are shown for values of 2 and 5, as discussed in the text. 

FIG. 2. Theoretical results for the [100] Compton profile in Li are compared for three 
different models: LDA without LP correction (solid), LDA with LP correction (light dashed), 
and the AGP scheme discussed in the text (heavy dashed) . 

FIG. 3. Effect of correlation on the Compton profile along two different directions in Be. 
The sohd curve gives the residual difference, AJ^'*^"^^ = Jagp — Jlda, between the LDA 
and the AGP profiles; note, that Jlda here includes the Lam-Platzman correction, which is 
also shown separately (dashed lines) for reference. 

FIG. 4. 3D momentum density along [110] is compared in Al for the LDA and the AGP 
case. 
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